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I. INTRODUCTION 

The inelastic (Raman) scattering of light by a semicon- 
ductor is an optical process that has proven its usefulness 



as a spectroscopic tool to investigate elementary excita- 
tions in semiconductors^^. From the theoretical point 
of view, the scattering process is described by a simple ex- 
pression coming from second-order perturbation theory 4 : 
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A fi is the quantum mechanical amplitude for the transi- 
tion from the initial (electronic) state, of energy Ei, 
to the final state, |/). This transition involves a change 
in the state of the radiation field. Indeed, the final state 
of the electron-photon system, \f,Ni — 1,1/), contains 
Ni — 1 incident photons of energy hv\ (one less than the 
initial state), and one photon of energy hvj (the scat- 
tered photon). The sum in Eq. JQ| runs over all inter- 
mediate (virtual) states. H e _ r is the electron-radiation 
interaction Hamiltonian, and Yi n t is a phenomenological 
damping parameter. 

From the amplitudes Afi, one computes the differential 
cross sectior 
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^AfifSiEi + hvi-Ef-hvf), (2) 



where dflf is the element of solid angle related to the 
wave vector of the scattered photon. Energy conservation 
is expressed by means of the delta function in Eq. J2J, 
which is approximated by a Lorentzian: 
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In the present paper, we focus on the Raman scattering 
in zero magnetic field from a quantum dot containing 
dozens of electrons. Thus, |z) and |/) are states with 
N e electrons. The incident laser energy, hvi, is taken to 
be resonant with an interband transition. It means that 
only the resonant contribution to A fi is considered in Eq. 
l|T)l£ and that the intermediate states, \int), contain an 
additional electron-hole pair. 
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Eqs. flj and J3J look very simple, but in fact their eval- 
uation is a cumbersome task because reliable approxima- 
tions to the many-particle wave functions \i), \int) and 
I/} need to be computed. A widely used simplified ex- 



pression is obtained by assuming a constant denominator 
in Eq. (|TJ and using completeness relations for the in- 
termediate and hole states. In this way, we arrive at the 
off-resonance approximation!: 
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where q and e* are the wave vector and the light polar- 
ization vector, respectively, a and a 1 label the Hartree- 
Fock (HF) states for electrons, and e and are elec- 
tron annihilation and creation operators. Notice that, in 
this approximation, the intermediate states play no role 
and the Raman amplitude is identified with the structure 
functions, i.e., only collective excitations in final states 
are supposed to contribute to the Raman peaks. Four 
terms are distinguished in Eq. The first one, pro- 

portional to e, ■ £j, corresponds to charge-density exci- 
tations (CDE). The next three, proportional to Ei x Sf, 
correspond to spin-density excitations (SDE). 

Most of the analysis of Raman experiments in quan- 
tum wells (qwells), wires (qwires) and dots (qdots) are 
based on expressions like Eq. (@J in spite of its known 
limitations. Experiments in qwells and qwires under ex- 
treme resonance (i.e., when the incident laser energy is 
close to the energy of the exciton) have revealed Raman 
peaks associated with single-particle excitations (SPE) 6 . 
These peaks do not arise from Eq. Q and are known 
to be related to taking a proper account of the inter- 
mediate (virtual) states!. For still higher excitation en- 
ergies (i.e., 40-50 meV above the band gap) a resonant 
enhancement of Raman intensities for particular values 
of hvi have been observed^. This effect is clearly not de- 
scribed by Eq. (0J . It has been ascribed to the existence 
of incoming and outgoing resonances in the intermediate 
states, although the nature of the outgoing resonances 
is not completely understood. The authors of Ref. Q 
have suggested the presence of higher-order Raman pro- 
cesses to explain the observed resonances. We shall show 
that the usual second-order expression Eq. (JTJ with a 
phenomenological Ti nt accounts for these effects. 

A review of relevant experimental facts of electronic 
Raman scattering in qdots can be found in Ref. |9J. In 
our opinion, the best experimental results are those re- 
ported in Ref. |lCj . As hvi moves from extreme reso- 
nance to 40 meV above it, the observed Raman spec- 
trum evolves from a SPE-dominated one to a spectrum 
dominated by collective excitations. The positions of col- 
lective excitations for the dots studied in Ref. [H| have 
been computed in Ref. |Tlj by means of expressions anal- 



ogous to Eq. but the dependence on hvi could only 
be obtained if one starts from Eq. QJ. 

In the present paper, we give a consistent theory of 
Raman scattering in medium-size qdots (dozens of elec- 
trons) based on the exact expression given in Eq. l(T]l. 
The theory is, however, "semiquantitative" because ran- 
dom phase approximation (RPA) - like wave functions 
and phenomenological Ti nt and Tf are used. The main 
limitation of the RPA functions in the present context is 
not related to the absence of correlation effects, but to an 
inadequate description of the density of energy levels and 
of the matrix elements of the electron-radiation interac- 
tion Hamiltonian. The main virtue of the RPA functions, 
on the other hand, is that collective excitations are de- 
scribed quite well. In spite of its limitations, the theory is 
able to reproduce all of the observed qualitative features 
of Raman scattering in qdots. 

With respect to previous calculations, we are aware 
of the exact computations (i.e., numerically exact elec- 
tronic wave functions plus Eq. Q) for a two-electron 
quantum ring made in Ref. |l2( , and of the approximate 
calculations for the 12-electron dot in Ref. In this 

paper, we report calculations for a dot with 42 electrons. 
Coulomb interactions are treated exactly (to the extent 
that the RPA approximation allows it) both in interme- 
diate and final states. Valence-band mixing effects for 
the hole are accounted for in the framework of the Kohn- 
Luttinger Hamiltonian. 

The plan of the paper is as follows. In the next section, 
we derive the theory needed for calculating the resonant 
Raman spectra, which are later discussed in section ITTT1 
Conclusions are drawn in section Hvi 



II. THE THEORY 

The computation of the Raman amplitude Afi re- 
quires: (a) the calculation of HF single-particle states for 
electrons and holes, (b) obtaining the final 7V e -electron 
states, I/), by means of the RPA scheme, (c) obtaining 
the (A^ e + l)-electron plus one hole states, \int), by means 
of the so-called particle-particle RPA formalism and, fi- 
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nally, (d) the computation of the matrix elements of the 
electron-radiation Hamiltonian, H e - r - Additionally, we 
shall compute matrix elements of multipole operators, 
something equivalent to the structure functions, and the 
density of final-state energy levels. Many of the required 
expressions and formulas were given explicitly in Ref. Q 
for the neutral electron-hole system. They can be used 
in the present context with minor modifications. 

Generally speaking, we use a HF-like scheme to de- 
scribe the ground state of the iV e -electron system (in fact, 
the RPA assumes that there are some "correlations" in 
the ground state, \i), as can be seen from the formulas 
below). An effective (conduction) mass approximation 
is used to describe electrons in the qdot. Thus, the N e - 
electron problem with confinement and Coulomb interac- 
tions is solved in this way. The excited states of this sys- 
tem, |/), are looked for with the help of the RPA ansatz, 
which has the form of a linear combination of "one parti- 
cle plus one (conduction band) hole" excitations over the 
ground state. To construct the intermediate states, \int), 
we need the HF (valence band) hole states. The latter 
are obtained by solving the Kohn-Luttinger Hamiltonian 
in the presence of the external confinement and the N e - 
electron background. The RPA ansatz for the interme- 
diate states has the form of a linear combination of "one 
electron (above the Fermi level) plus one (valence band) 
hole" excitations over the ground state. As a direct result 
of the RPA calculations, we obtain matrix elements such 



as (f\e^.ix\i) , which are needed to compute the Raman 
amplitude. 



A. HF states for electrons 

We will model the qdot with a disk of thickness L. 
The disk axis coincides with the z axis. At z — and 
z = L a hard wall potential confines the electron mo- 
tion. On the other hand, the in-plane confining potential 
will be assumed to be parabolioi^, with a characteristic 
energy fno^. The HF electron single-particle states (or- 
bitals) are expanded in terms of two-dimensional (2D) 
oscillator wave functions^, ip S) according to the follow- 
ing ansatz: 




(k?irz/L)J2cr<Ps(x,y)Xs(<r), (5) 



where L is given in nanometers, k z is the sub-band label, 
and Xs are spin functions. 

The expansion coefficients, Cf , and the energy eigen- 
values, E a , are obtained from a set of equations similar 
to Eq. (6) of Ref. [l4|, in which hole contributions shall 
be ignored: 
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where 7 runs over occupied states (/xp is the electron The HF equations are solved iteratively. Twenty oscil- 
Fermi level), and the 2D oscillator energies, in meV, are lator shells are used in the calculations, 
written as: 



To be definite, we will use parameters appropriate for 
GaAs, i.e., the conduction band effective mass for elec- 
trons is m e /mo = 0.067, and the relative dielectric con- 
stant is k — 12.5. 

Two-dimensional Coulomb matrix elements^ will be 
used instead of the truly 3D ones. Consequently, we will 
assume that the matrix elements will be diagonal in the 
sub-band index, k z , and will multiply the matrix ele- 
ment by a strength coefficient, f3, in order to simulate 
the smearing effect of the z directional. The coefficient 
will take values from 0.6 for L = 25 nm, to 0.8 for L = 8 
nm. 



B. HF states for holes 

To guarantee that both electrons and holes are con- 
fined in the same spatial region, we will assume different 
confining potentials, i.e., we will require: 

m e uo = m\ h uo^ h = mfuj l Q h , (8) 

where m^ h is the in-plane mass of the j = 3/2, rrij — 
±3/2 (heavy) hole, and m 1 ^ is the mass of the j = 3/2, 
nij = ±1/2 (light) hole. The ansatz for the HF hole 
orbitals is the following: 
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The expansion coefficients and energy eigenvalues are to 
be determined from the equations: 

= Vz £ sin ) 

s,k Z) mj 

(9) 
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The first term is the Kohn-Luttinger Hamiltonian, 
"Hjcc, whose matrix elements are given in Appendix lAl 
The second term is the electrostatic field of the back- 
ground N e electrons. Coulomb interactions are assumed 
to be diagonal in rrij indexes as well. Notice that, be- 
cause of the form of TLkc , hole states are grouped into 
sets with a common value of fh = —"m^ + 1st where / s is 
the angular momentum quantum number corresponding 
to the hole oscillator state (p s . 

HF energies and wavefunctions for electrons and holes 
are used as input in the calculation of the many particle 
functions |/) and \int). 

C. The final states 

The final states, | /) , are excitations of the iV e -electron 
system. In the RPA, they are obtained as linear combi- 
nations of "one particle plus one hole" excitations over 
the initial state, 



l/}= E (XrxeUx-Y^eie^li), (11) 

where the index A runs over occupied HF states, and 
a runs over unoccupied states. Detailed equations for 
the coefficients X, Y and the energy eigenvalues can be 
straightforwardly obtained from the formulas (12-15) of 
Ref. ]l|: 

/J {AtA.tm^tm + BcrX^rYfir} = ^fX^x, 
/x H~ -1 Act, /at 

in which hflj is the excitation energy, r and fj, are in- 
dexes similar to cr and A, respectively, and the A and B 
matrices are given by: 
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Notice that X^x has a straightforward interpretation 
in terms of a transition amplitude: 

X IT x = {i\e t xe IT \f), (14) 

and similarly for the Yxrr- 

We shall stress that final states are characterized by 
the quantum numbers Al and AS Z representing the vari- 



ation, with respect to the ground state, of the total angu- 
lar momentum projection and the total spin projection, 
respectively. Conventionally, we will call Al = states 
as monopole excitations, Al — ±1 states as dipole exci- 
tations, Al = ±2 as quadrupole excitations, etc. 

The calculation of strength functions defined in Eq . 
follows also very simply from the results of Ref. [14( . 
One expands the exponential 
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e i{qi-qf)-f _ e i(q*i-qzf)z ^(q\\i-q\\s)-P 

(15) 

and makes use of the definition of multipole operators, 
d l given in Ref. Q: 



(a|p 2 | 7 ); 1 = 0, 



(16) 



whose detailed expressions can be found in Appendix B 
of that reference. In the later formulas, a denotes the 
orbital part (no spin function included) of the HF elec- 
tronic state a. The spin projection quantum number is 
explicitly indicated in Eq. (@J. With respect to the part 
depending on z, one uses that 



(k z \e lqz \k' z ) 



Aik z k' z q (-1 + e i7r « cos(fc z 7r) cos(fc») 
tt [{k z - k' z y - f] [(k z + k' z y - 



(17) 

The strength functions, or more precisely the multi- 
pole operators, allow a further classification of the |/) 
states into collective and single-particle excitations^. A 
charge monopolar collective state |/), for example, gives 
a significantly nonzero value for the matrix element: 



\int)= J2 V^elhUi), (20) 

<t>Mf,t 

where cr is a HF electron state above the Fermi level, 
and r a HF hole state. The equations for the expansion 
coefficients, V , and the energy eigenvalues are explicitly 
written in Ref. 
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The quantity Ml int gives the excitation energy, measured 
with respect to the ground state of the A^-electron sys- 
tem, and the coefficients Vo~r can be interpreted as the 
transition amplitudes: 



V aT = {i\h T e„\int) . 



(22) 

The intermediate states are characterized by the quan- 
tum numbers 



•F — J hi S z , 



(23) 



where l e and S z are the angular momentum and spin 
projection of the added electron. 



E. The geometry of the Raman experiment 



(18) 



whereas for a single-particle excitation, the matrix ele- 
ment practically equals zero. By "significantly nonzero 
value" we mean that |£>^J 2 is greater than 5% of the 
energy- weighted sum rule for the monopole operato»i2ii&: 



f 



2h 2 
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Similar criteria can be formulated for charge multipo- 
lar states^, or for spin excited states (involving or not 
spin reversal with respect to the ground state). The lat- 
ter are related to the last three terms of Eq. J2J. 



D. The intermediate states 

The intermediate states with N e + 1 electrons and one 
hole, can be obtained from the so-called particle-particle 
Tamm-Dankoff approximation (pp-TDA), which is an 
uncorrelated pp-RPA function, i.e., no particles below 
the Fermi level are created^: 



In the present paper, we restrict ourselves to the so- 
called backscattering geometry, which is often used in 
experiments. The incident laser beam is deflected inside 
the dot because of SnelPs law. Thus, the actual angle of 
incidence (with respect to the z axis) is: 



arcsm 



sin < 



V 



(24) 



where r\ « 3.5 is the qdot refractive index. If qi denotes 
the wave vector in vacuum, then inside the dot we have: 



<1\\ 

and for the scattered light: 



qi sin < 



l'\\f 



q f sin (j) f , 



(25) 
(26) 



(27) 
(28) 



-r}q f cos (p f , 
with (pi = (pf. 

We distinguish between the polarized geometry, in 
which e*j and £/ are parallel (along the y axis in our cal- 
culations) , and the depolarized geometry, in which e f (in 
the xz plane) is orthogonal to Ei. A detailed definition of 
angles for more general geometries can be found in [T^j . 
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F. The matrix elements of H e -r 

ii~_ r is the part of the electron-radiation Hamiltonian 
responsible for the annihilation of an incident photon and 
creation of an e-h pair. Its matrix elements are written 
asii: 

(int\H-_ r \i) ~ (band- orbital)®. V* T , (29) 

tr>fJ, F ,t 

I 



The band factor, Si ■ po- 3 , mj , is computed according to 
Tabled The factor (k z \e lq " z \k' z ) is computed with the 
help of Eq. H17|) . Finally, the computation of the integral 
(the orbital factor) is made along the lines sketched in 
Ref. Q. 

On the other hand, H^-r * s that part of the electron- 
radiation Hamiltonian responsible for the creation of a 
photon and annihilation of an e-h pair. Its matrix ele- 
ments are given as 1 *: 

(f\H+_ r \int) ~ ]T J2( band -° rbltal ^V^Kx, 

(31) 

where Xa-x is one of the coefficients entering the RPA 
expansion Eq. Qllfl. The band-orbital factor is obtained 
from Eq. (|30[) by replacing i by / and taking the complex 
conjugate of the expression. 

G. Phenomenological Tf and Tint 

The main decay mechanism of electronic excited lev- 
els in a qdot at very low temperatures is the emission 
of longitudinal optical (LO) phononsi^. We will ignore 
surface effects in a qdot, and assume a threshold excita- 
tion energy appropriate for GaAs, Tholo ~ 30meV, for 
the emission of LO phonons. 

Only final states with excitation energies lower than 
Swlo will be considered in order to exclude Raman peaks 
related to phonon excitations. It means that final states 
will have small widths, for which we will take a constant 
value, Tf, in the interval between 0.1 and 0.5 meV. 

In the same way, for intermediate states with exci- 
tation energy lower than Hlolo we will take Tint =0.5 
meV. For higher excitation energies, the LO phonon de- 
cay mechanism becomes active and the widths suddenly 
increase. In this case, we will take Ti nt =10 meV, except 
for a set of particular states, which can be interpreted 
as "excitons" or "excitons + plasmons" , whose meaning 



where Vo-t are the coefficients entering the pp-TDA Eq. 
(|2L)p. Because of valence-band mixing, the orbital (en- 
velope) and band wave functions of the hole get mixed. 
The band-orbital factor in Eq. I|29|) is defined as: 



(30) 

I 

will become clear below in the discussion of Raman scat- 
tering well above band gap. In this latter situation, we 
will take Ti nt —2 meV. 

We stress that the role of Tf and T^t as functions of 
the excitation energy in the Raman spectra has not been 
pointed out before. In our view, the qualitative change 
of the resonant Raman spectrum when the incident laser 
energy, hvi, is raised 30 or more meV above the band gap 
is related to the sudden increase of Ti n t . 



III. CALCULATED RAMAN SPECTRA 

In the following, we report results for a 42-electron dot. 
The disk thickness and the harmonic confinement take 
values L = 25 nm, and Tiloq — 12 meVpS respectively. 
The chosen Hcoq corresponds to a qdot in the strong con- 
finement regime, and the number of electrons to a closed- 
shell quantum dot. 

We show in Fig. ^the electronic excitations of the dot, 
i.e., the spectrum of final states, |/). The reference en- 
ergy is the energy of the ground state, \i). The excitation 
energy is precisely what is measured as the Raman shift 
in the experiments. 

To the left of the y axis, states with AS Z = 1 (with 
respect to \i)) are represented, while to the right of the 
y axis, states with no spin flips are shown. In the fig- 
ure, we identify the collective excitations, labelled CDE 
and SDE, and give explicitly the corresponding fraction 
of the energy-weighted sum rule 14 . In the Al = 0, 
AS Z = case, for example, the CDE state concentrates 
the strength of the charge monopolar transition (from |i) 
to I/)), and the SDE state concentrates the strength of 
the spin monopolar transition (with no spin flip). The 
rest of the states shown correspond to single-particle ex- 
citations (SPE). Notice that, in general, collective exci- 
tations are isolated from SPEs. 

For the intermediate states, we take a nominal band 
gap, Eg Cl p = 1560 meV. This gap is renormalized by 
Coulomb interactions. We define the renormalized E' iri 



(band- orbital)® =J2 £ C£gX (ej -jv^.) (k^^K) j ' e^MMMp) > 

s,k z t,k' z ,mj 
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as the energy of the lowest intermediate state. Note that 
this convention may not coincide with the experimental 
definition of the effective gap in terms of the position of 
the exciton line. 



A. Raman spectra below the effective band gap 

Measurements of electronic Raman spectra when the 
laser excitation energy, hfi, is below E' gap have not, to 
the best of our knowledge, been reported for qdots. In 
the present section, however, we show that such mea- 
surements could provide information for both collective 
excitations and SPEs in qdots. Raman intensities for 
both kinds of excitations show comparable magnitudes. 

Note that we use only the resonant contribution to Afi, 
Eq. , in spite of the fact that the present situation does 
not correspond, strictly speaking, to a resonant process^ 

We have the possibility of computing the spectrum for 
each multipolarity of final states. Results will be pre- 
sented in this way, although in an experiment all the 
multipolarities can be observed in the same spectrum. 

We show in Fig. Elthe polarized Raman spectrum for 
monopole final state excitations, computed with = 0.5 
meV. E'„ in this situation is 1599.2 meV. The incident 
(and backscattered) angle is equal to 20°. hi>i is swept in 
a 30 meV interval below E'„. Notice the monotonic in- 
crease of intensities as hvi rises. One peak corresponding 
to the CDE, and a second one related to the SPEs are ob- 
served. In the latter case, there is a group of energy levels 
contributing to the peak in the figure. We may think of 
this set of levels as a Coulomb-renormalized oscillator 
shell. As we are dealing with monopole excitations, the 
average position of the group should correspond to 2hio' 0l 
i.e., the renormalized oscillator energy is huj' « 9 meV. 

The fine structure of the SPE peak is shown in Fig. 
along with the density of final-state energy levels. In 
this case, the monopolar polarized Raman spectrum was 
calculated with T f = 0.1 meV. The depolarized spectrum 
is shown in Fig. [3p (In fact, only the energy interval 
corresponding to the SPE peak is shown. The SDE peak 
is outside this interval). Histograms with a step of 0.1 
meV are used to represent the level density. Although 
Eq. Q refers to collective excitations, we have used its 
implications to correlate the polarized Raman spectrum 
with the level density of charge monopolar SPEs, and 
the depolarized spectrum with the density of monopolar 
spin excitations. The Raman spectra reproduce quite 
accurately the details of the level density in both cases. 
A general remark concerning Fig. [3|is that the intensity 
of the depolarized peaks is about three times lower than 
the intensity of the polarized ones. 

Depolarized spectra for spin-flipped monopolar and 
dipolar final states are shown in Figs. 0^ and b. Polar- 
ized spectra for dipolar and quadrupolar states are shown 
in Figs. 0J: and d. The shell structur o 10 i 21 is clearly seen 
in these figures. Dipolar and spin-flipped final states are 
strongly depressed in the Raman spectra. Quadrupole fi- 



nal states, on the contrary, show magnitudes comparable 
to monopolar peaks. 

The fact that even multipoles are favored in the Raman 
spectra is understood in terms of the even parity of final 
states in a two-photon process. On the other hand, spin- 
flipped final states are reached only as a consequence of 
valence-band mixing. We (virtually) create an electron 
and a hole, whose dominant component is rrij, and anni- 
hilate the same hole and an electron with opposite spin. 
The amplitude of the latter process is proportional to the 
minority component of the hole wave function, Xm'- ■ It 
means that the amplitude squared |j4/i| 2 will be propor- 
tional to \x m r\ 2 - 

These calculations show that experimental measure- 
ments of the electronic Raman spectra with below band 
gap excitation can provide valuable information on the 
collective states and SPEs of qdots. Furthermore, below 
band gap excitation can overcome the problem of overlap 
with the intense photoluminescencc observed under reso- 
nant excitation. The peak maxima exhibit a continuous 
but not very marked increase in intensity with excitation 
approaching the band gap (see Fig. |5J, indicating that 
excitation around 30 meV below the gap is sufficient. The 
other notable feature of these calculations, apart from the 
marked differences predicted in Raman intensities of the 
polarized and depolarized multipolar components, is the 
fine structure of the SPE Raman peak (see Fig. |3J). It 
would be interesting to probe all these aspects experi- 
mentally. 

B. The extreme resonance region 

In the present section, we consider Raman scattering 
when hvi moves in a 30 meV window above E' gap . We 
will call this interval the "extreme resonance" window. 

Fig. [SJshows a polarized Raman spectrum correspond- 
ing to charge monopolar final states. As in Fig. [2 we 
used T f = 0.5 meV. Two characteristics of Fig. [21 make 
it very different from Fig. [21 (i) Peak intensities are not 
monotonous with respect to variations of hvi, and (ii) 
The position of the maximum in the SPE peak moves 
slightly with hvi. Both properties are related to reso- 
nances in the intermediate states. 

Resonances in the intermediate states can be better vi- 
sualized if we follow the Raman intensities of the peaks 
shown in Fig. [3J For this purpose, we computed monopo- 
lar Raman spectra with Tf = 0.1 meV and varying hvi 
with a 0.5 meV step. The results are drawn in Fig. 
The monotonous increase of peaks in the hvi < E' gap re- 
gion is apparent in the figure. On the other hand, for 
hvi > E' gap the intensity variation with laser excitation 
energy is much more complicated. The intensities of the 
individual SPE components rise and fall markedly with 
laser energy, as has been observed experimentally. This 
variation is attributed to individual resonances occurring 
within intermediate states lying close to the band gap as 
the incident light energy sweeps through them. The ef- 
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feet is particularly noticeable for hvi ~ 1616 and 1626 
meV. In the associated intermediate states, the added 
e-h pair has zero total angular momentum projection, 
and the hole is basically a heavy hole. Notice that the 
same intermediate states are responsible for the strong 
enhancement of Raman intensities in both the polarized 
and depolarized geometries. 

A spin monopolar SPE Raman peak is followed as a 
function of hvi in Fig. Qd. For comparison, we have 
also given the products \(f\H^_ r \int)(int\H^_ r \i)\ 2 for 
each intermediate state, and the absorption strengths 
\(int\H~_ r \i)\ 2 (the upper panel). The optical absorp- 
tion coefficient is defined according to 

« = ElMt,Bf |e .^ tI .^ (32) 

Fig. shows that, in the present situation, peaks in the 
optical absorption coincide with peaks in the Raman in- 
tensities, which leads to the conclusion that the latter are 
related to incoming (absorption) resonances. Note that 
from Fig. [7jp it follows that interference effects among 
intermediate states are weak in the Raman scattering un- 
der extreme resonance^. Raman spectra in the extreme 
resonance region look similar to the spectra shown in Fig. 
0] but with SPE peaks much higher than collective ones 
and selectively enhanced for particular values of hvi. A 
comparison with the density of energy levels would lead 
to results very similar to those of Fig. |3J 



C. Raman spectra 40 meV above band gap 

As mentioned in Section HUG!, we assume that T int ex- 
periences a sudden increase when Ei nt > E' gap + Hujlo- 
As a result, the contribution of these states to A/<, Eq. 
JTJ, loses its resonant character even when hvi sweeps 
this energy range. It means that peak intensities become 
smooth functions of hvi, as for below band gap excita- 
tion. Both collective and SPE Raman peaks decrease in 
intensity for hfi > E' gap + hwho ( a s compared with val- 
ues at extreme resonance), but the SPE peaks are more 
strongly depressed^. Fig. |S]shows a typical spectrum at 
hvi = 1642 meV. 

Nevertheless, a modest increase of the peak intensity 
for hvi well above the band gap may result not only from 
large values of \(f\H^_ r \int) (int\H~_ r \i}\ 2 but also from 
relatively small level broadening (as compared with the 
neighboring levels). There could be a set of intermediate 
states in which Ti nt takes relatively small values. One 
can think, for example, of the lowest state in a sub-band, 
let us say the k z = 2 electron sub-band. Inter sub-band 
transitions due to phonon emission are not as fast as 
intra sub-band transitions 19 . Consequently, the Ti nt of 
the lowest state in the sub-band is relatively small. We 
will call these states "excitonic" states, X. If the product 
(f\H^_ r \X)(X\H~_ r \i) is not small, a peak in the Raman 



intensity will appear. In the interpretation of Ref. 0, this 
is an "incoming" resonance. 

On the other hand, "outgoing" resonances correspond 
to emitted photons with energy Ex, as shown in Fig. [5] 
The intermediate states are located at excitation energies 
Ex + Ef. One can think of these states as an exciton on 
top of a collective excitation or, conversely, a collective 
electronic excitation on top of an excitonS. The X in 
this case may correspond to an absorption peak in the 
extreme resonance region or well above band gap. 

The resonant enhancement of these "exciton plus plas- 
mon" states can be due, again, to big numerators or small 
denominators in Eq. (JTJ. Relatively small Ti nt could be 
related to the collective nature of these states or to the 
relative isolation from neighboring levels. 

To illustrate the effect of resonant enhancement of Ra- 
man peaks, we show in Fig. ^|the intensity of the CDE 
peak as a function of hvi. The entire range of variation is 
shown for completeness, i.e., below band gap excitation, 
extreme resonance, and well above band gap excitation. 
In the extreme resonance case, the enhancement is re- 
lated to absorption maxima as mentioned above. On the 
other hand, for hvi well above band gap, we pick up an 
intermediate state with energy Ex + Ecde, where the 
exciton level X corresponds to the absorption maximum 
at 1620 meV (see Fig. GJi), and E C de ~ 22 meV is the 
energy of the charge monopolar collective state. For this 
state, we chose Ti nt = 2 meV. The effect is an enhance- 
ment of the intensity around 1642 meV, as shown in the 
figure, that according to Ref. |8jis an outgoing resonance. 

This calculation reveals that the dependence of Ti nt on 
Ei n t may dictate the qualitative features of Raman scat- 
tering with well above band gap excitation. A consistent 
treatment in which the Ti nt are computed for each \int) 
is left for future study. 



IV. CONCLUSIONS 

This theoretical investigation of the role of resonance 
excitation in electronic Raman scattering from qdots has 
revealed many hitherto unsuspected features and details. 
In general terms, the Raman intensities of the SPEs, 
CDEs, and SDEs are strongly affected when the incident 
light energy is swept from below, through, and above 
the quantum dot band gap. Incoming resonances pro- 
duce a rapid variation in intensity for excitation energies 
just above the band gap, and outgoing resonances are 
predicted for higher excitation energies, as observed ex- 
perimentally. In fact, observation of the Raman intensity 
of just one SPE, for example, as a function of the inci- 
dent light energy can provide precise details of the optical 
absorption spectrum and density of states. The role of 
damping in the intermediate states has been shown to be 
a significant factor in determining these resonances and 
deserves further theoretical analysis. 

Another aspect of this work is the unravelling of the 
complexity of features in polarized Raman spectroscopy 
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of qdots. This spectral complexity in dots with large 
numbers of electrons has been evident from the first ex- 
periments. These calculations have shown what excita- 
tions dominate in which polarizations and point the way 
to a better control of what is measured in future experi- 
ments. 

Experimentally, the predicted selective resonances 
have advantages and disadvantages. By employing reso- 
nance excitation, a particular final state can be enhanced 
over its companions and thus make it easier to iden- 
tify. On the other hand, the rapid variation in inten- 
sity of the individual and numerous SPEs makes it dif- 
ficult to uniquely identify them. In this regard, the fine 
structure evident for SPEs from these calculations should 
be explored by performing high resolution Raman spec- 
troscopy in future. The best situation for evaluating the 
SPEs would be for an excitation energy about 30 meV 
below the band gap, where some overall resonance en- 
hancement occurs but the contribution from band gap 
photoluminescence would be weak. No such experiments 
have been performed to date. 
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APPENDIX A: MATRIX ELEMENTS OF THE 
KOHN-LUTTINGER HAMILTONIAN 

In the present appendix, we give the matrix elements 
of the Kohn-Luttinger Hamiltonian entering Eq. i|l(J|l . 
Results are presented for the more general case when a 
magnetic field, B, is applied in the z direction. We use 
the following set of parameters^: 



71 = 6.790, 72 = 1.924, 73 = 2.681, 
7=(72 +73)/2, k = 1.2, g = 0.04 

The Kohn-Luttinger Hamiltonian has the following 
structure in the m,j variable: 



( H 3 i 2 S 

S { H 1/2 
fit 

V 



R 
R 

H_ 1/2 -S 
i?t -St H _ 3/2 J 



(Al) 



The H terms are diagonal in oscillator and k z indexes. 
They are given by: 



H±z/ 2 = Ml, 



TO e (7i +72) 



m 



(2fc+|I h | + l) + 



h 2 , kW 
2^~ 2 ^ 



fiWce (7l +72) 



mo 



lh ± 3k 



H 



±1/2 



; (7i " 72) 
too 



(2fc+|/ fc | + l) + — (71 + 272) 

ZTOo 



k 2 z ir 2 hu; ce m e {^i - 72) 



L 2 



mo 



lh ± V-B 



(A3) 



where too is the electron mass in vacuum, /ib is the 
atomic Bohr magneton, cu ce is the electron cyclotron fre- 
quency, f2 e = y/uj 2 + ^ e /4, and k, lh are the radial and 
angular momentum quantum numbers corresponding to 



the 2D oscillator state. 

The matrix elements of the S and R operators are writ- 
ten in the following form: 



(k,lh,k z \S\k ,l h ,k z ) — 



\/37i 2 73 



too 



h L 



7T3 7757 $1 h . I' + 1 &k z + fet , odd 

{k 2 -k 2 ) h h 



VF+T (1 - w ce (2ft e )) 4,fc<+i + y/kf + \l' h \ (l + w ce /(20 e )) 4,fc' ; l' h <-l 



(A4) 



/ F(l+Wce/(2ne))*fc,fc'-i + -«W(2fi e ))4,fc'; l 'h > 0, 



(k,l h ,k z \R\k',l' h ,k' z ) = —^-—^hn e 5 lh ^ h+2 5 kz . K 
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x < 



' v/( fc ' + Kl - !)( fc ' + K\) C 1 + u ca /(2n e )) 2 5 Kk , 

2y/(k' + l)(k' + \l' h \) (l - ( Wce /(2fi e )) 2 ) 4,fc'+i; l' h < -2 



x < 



[ y/(k' + 2){k' + 1) (1 - u> ce /(2Q e )y S k , k . 

' -y/k'(k' + 1) (1 + LO ce /(2n e )f 

-2fc' (l - (^ ce /(2f! e )) 2 ) S kM ; 



+ 2 



^ = "I 



{ -y/tf + 2){k' + 1) (1 - w ce /(2fi e ))^ 5 k , k 



'+1 



x < 



' v /(fc'-i)fc'(i + ^ ce /(2r! e )) 2 4^-2 

2 v /fc'(fc' + + 1) (l - (co ce /(2ft e )) 2 ) 4,fe'-i; ^>o 
[ + l' h + 2)(k' + V h + 1) (1 - Lu ce /{2Q e )f 6 k , k >. 



(A5) 
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TABLE I: The quotient e, ■ p a , mj /{iP), where P is the GaAs band constant. 
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FIG. 1: Spectrum of final states in the model qdot. 
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FIG. 2: (Color online) Polarized monopolar Raman spectrum below band gap (Tf = 0.5 meV). 
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FIG. 3: (Color online) Polarized and depolarized monopolar Raman spectra (Tf = 0.1 meV) and comparison with the density 
of energy levels. 
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FIG. 4: Raman spectra in different channels (Ff =0.1 meV). The incident laser energy is 1594.2 meV. 
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FIG. 5: (Color online) Monopolar polarized Raman spectra (Tf = 0.5 meV) in the extreme resonance region: E' gap < hvi < 
E' + 30 meV. 
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FIG. 6: (Color online) The dependence on hvi of the intensity of SPE peaks identified in Fig. 
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FIG. 7: (Color online) (a) Absorption in the model qdot. (b) Intensity of the spin monopolar Raman peak with AEf = 18.6 
meV as a function of hvi. 
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FIG. 8: Raman spectrum in the polarized geometry (Tf = 0.1 meV) for well above band gap excitation. 
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FIG. 9: Schematic representation of incoming and outgoing Raman resonances. 
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